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Introduction. The viscoelastic behavior and viscous-flow properties of in- 
compressible materials undergoing flnite deformations can be simulated in terms 
of models of integral and differential type. The integral models describe the his- 
tory of medium deformation using integral equations. However, these equations 
are difficult to generalize to complicated problems such as, for instance, simula- 
tion of mass exchange processes in dissipative media. Attempts to present dam- 
age accumulation and thixotropy of elastomeric materials in terms of integral 
models add complexity to these models and hamper their practical application. 

The differential models describe the rheological properties of materials in 
terms of tensor internal variables and, as a rule, assign them the physical mean- 
ing of stresses PEHSHiEllSlEI or strains IBIEIII01IIIIII21II31II1I • These models can 
be conveniently represented with the aid of symbolic diagrams illustrating the 
mechanical behavior of the medium. In the paper presented here a fairly simple 
method of obtaining constitutive equations is described. The method has been 
constructed using symbolic diagrams. The main advantage of the method is 
rather simple physical meaning of mathematical expressions. The method has 
its origin in the ideas reviewed in works flSlIISlIIZI • 

Singular points in the symbolic diagram. This work presents the 
method allowing one to construct the mathematical model for the medium 
whose mechanical behavior is illustrated by a diagram consisting of blocks of in- 
terrelated elastic and viscous elements. Consider this method using the diagram 
shown in Fig. [T] as an example. To each junction of elements from the symbolic 
diagram illustrating the material behavior, a vector in nine-dimensional vec- 
tor space corresponds. The space forms a set of deformation gradient tensors. 
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Fig. 1: Schematic representation of the model for viscoelastic material Arrows show hnear 
mappings of the infinitesimal domain of the material. Corresponding deformation gradients 
present information about the mappings which couple reference and examined points in the 
diagram and are used for introducing internal variables (point coordinates of the diagram in 
nine- dimensional vector space) 



Between the points there are elastic and viscous interrelations. Special role is 
assigned to three points: left, right and reference. 

The left point in the scheme is denoted by symbol A (Fig. [T]) to which the 
deformation gradient of the material corresponds. That is, 



Fa = F = 



dx 
dxo 



where x = x{t, xq) is the radius-vector of location of medium points at current 
time t; Xq is the radius-vector of location of medium points at reference time. 
It is evident that the rate of deformation tensor D and the spin employed 
in the theory of nonlinear media, are connected with point A. 

The right point in the scheme is denoted by symbol £ in Fig.[Tl Its peculiarity 
is that the strain gradient Fg, corresponding to point is the rotation tensor 
which determines the rotation of the infinitesimal bulk of the medium as 
an absolutely solid body and is independent of the material deformation 



Fs — Rm- 



A 



Fig. 2: Diagram of the ith. symbolic element. In constructing of constitutive equations, defor- 
mation gradients i^^, Fq and tensor measures of the rate of deformation of Z)^, Z)g 
are assigned to points A and B. 

Hence, the rate of deformation tensor, calculated for the right points of the 
system, is equal to zero 



The reference point is denoted by symbol J- in Fig. [H In nine-dimensional 
space, the unit tensor / corresponds to this point 



Tensors used for establishing the strained state of the elements 
in the symbolic diagram and the rate of change of this state. Let 

us examine some peculiarities of the mathematical description of the strained 
state of the arbitrarily selected element of the system (Fig. EI). Assume that the 
element is numbered by and its strained state is determined by comparing 
two deformation gradients and two tensors of rate of deformation. 

The element i is connected with other elements at points A and B. De- 
formation gradients and correspond to points A and B. Using these 
gradients, the rate of deformation tensors D^, Db and spins VF^, Wb are 
determined as 




Ft 



= I. 







4 



To interpret tensor values from physical standpoint, it is reasonable to elab- 
orate the appropriate mathematical model of the medium starting from the 
following premises. The rate of deformation tensor of the ith element Di can 
be calculated for viscous and elastic elements from the formula 

D, = Da- Db. (1) 

Furthermore, in the elastic elements of the scheme the relation between defor- 
mation gradients, corresponding to points A and can be defined by the law 
of multiplicative decomposition 



A 



V, Fs, (2) 



where the stretch tensor of the ith element Vi is found by three orthonormalized 
eigenvectors n^^ n\ (defining the direction of principal tensor axes) and by 
stretch ratios A^, A2, A3 

3 

k=i 

It is easy to check that condition determines the relation between tensors 
Da, Wa and D^, Wb in elastic elements 

Db = I V;' B, {Da + 1^7) + {Da + Wa) B~ B, j Vt\ (3) 

-B,{Da + Wa^)^ V.\ (4) 

and the rate of time variation of the stretch ratio 

Al = Ai,n*fc®n^ A, (5) 
where the tensor Bi is the left Cauchy-Green tensor of the elastic element 

B, = V?. 

Coupling of elastic elements in the diagram. In the proposed method 
of constructing of constitutive equations there is a condition which should be 







fulfilled when coupling the elements. The elastic elements must be connected by 
their left ends with other elastic elements or with the left point in the symbolic 
diagram. This restriction is caused by the necessity to calculate by formula (01) 
the time derivatives of tensors Bi for all elastic elements using for this tensors 
Di, D and W. Calculation starts with elastic elements connected by their left 
ends with the left-hand point in the diagram. At this step we find not only 
the time derivatives of tensors Sj, but also the rate of deformation and spin 
tensors of the right junction with the adjacent elements (tensors Db and Wb in 
Fig. 121). Calculations are performed using formulas dH), (01), OH). Then we carry 
out calculations for elastic elements connected on the right with the elements 
just examined, and so on. 

Simulation of mechanical properties of medium elements. In isother- 
mal case, as an expression for mass free-energy density /, we use the function 
of stretch ratios A*^ of the elastic elements of the medium 

where n is the number of elastic elements. It is assumed that each elastic 
element satisfies the volume invariance condition when subjected to deformation 

AlA^A^ = 1. (7) 

Therefore, the Cauchy stress tensor is calculated by the formula from the non- 
linear elastisity theory for incompressible materials 

T,=pa+ pY,K§^ril0nl (8) 

where p is the medium density, pi is an indefinite parameter, and i is the 
elastic element number. The characteristic features of viscous elements can be 
simulated by ordinary state equations for viscous liquids 

Tj=p,I + 2rjjDj (9) 

and their incompressibility condition 

I Dj = 0, (10) 

where Tj is the Cauchy tensor for viscous elements, r]j is the coefficient of shear 
viscosity, and j is the viscous element number. 

Interrelations between stresses for coupled elements. For determina- 
tion of the closed system of constitutive equations, it is necessary to establish 
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the interrelations for stresses at diagram junctions. They must be formulated 
in such a way that the sum of all forces acting on the arbitrarily selected junc- 
tion would be equal to zero. Since we study the points in nine-dimensional 
vector space (space generates a set of possible values of deformation gradient 
tensors), then we can use the Cauchy stress tensors of the diagram elements as 
the nine-dimensional force vectors. 

Now let us formulate the rule. The sum of Cauchy stress tensors of the 
elements coupled on the left is equal to the sum of Cauchy stress tensors of 
the elements coupled on the right at any diagram junction. It is assumed in 
this case that the forces, determined in nine-dimensional vector space by the 
Cauchy stress tensor, act from left on the left point of the symbolic diagram 
and from right on its right point. It means that for the diagram given in Fig. [T] 
the following equalities are valid: 

T = Ti + T2 + T3 + T4 (at point A), 

Ti = T5 (at point B), 

T2 + T5 = Tq (at point C), 

T3 + T6 = T7 (at point D). 

The right point allows us to formulate the linearly dependent equation 

T4 + T7 = T (at point S) 

which gives no additional information. 

Dissipation inequality. In the isothermal case, the thermodynamic con- 
straint (second law of thermodynamics) imposed on processes in the medium 

T-D- p f>0 

takes the form 

n n+m 

1=1 j=n+l 

where elastic and viscous elements are numbered by i and j, respectively, n is 
the number of elastic elements in the system, and m is the number of viscous 
elements in the system. Finally, using the type of free energy (Ej) and equlities 
®' ©"(lini), the dissipation inequality can be written as 

n+m 
j=n+l 
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Fig. 3: Combination of elastic and viscous elements (a) allows description of large hysteresis 
losses at large deformations and slightly changes the material behavior at small deformations 
(b) 



This inequality is satisfied automatically provided that the shear viscosity co- 
efficients can not be negative. 

Example. Consider the material whose behavior is illustrated by a diagram 
in Fig. El Dots denote the blocks of elements. We will not describe here the 
way of coupling elements into these blocks. 

In modeling the viscoelastic properties of polymers, a proper consideration 
of hysteresis losses is of particular importance. Usually the most challenging 
task is to provide the accurate description of large losses in the region of large 
values of stretch ratios and small losses in the region of small deformations of 
the material. In the considered example it is possible to describe such material 
behavior using the potential of free energy 

./ ^ - E ^ in (l - ^ ) . . (/.(BO - 3) - 3) , 

where ci, C2, C3, J^, are material constants, and /i(-Bj) is the first invariant 
of tensor Bi. The first two terms in the potential are known in the literature 
as Gent potentials formulated for corresponding elastic elements flBj. 

Two peculiarities of the material behavior defined by elements under con- 
sideration should be noted. During the first loading cycle there are hysteresis 
significant losses and during the second and other cycles the losses are essen- 
tially less. However, such behavior of the material is viscoelastic by nature, and 
not a manifestation of MuUins softening effect. 



Another distinguishing feature of the model is that the free energy poten- 
tial may involve the cross terms interrelating the behavior of different elastic 
elements. In the example presented, this is the last term. It means that the 
deformation of the third element changes the elastic behavior of the first, and 
the deformation of the first element specifies the properties of the third. The 
physical meaning of symbolic elements remains unchanged. 
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